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Derivation of equations governing the smoothing, of 
rough surfaces, based on Mullins' (1957,1960,1963) 
theories of thermal grooving and of capillarity- 
governed solid surface morphology. As an example, 
the smoothing of a one-dimensional sine-shaped. surface 
is discussed. 
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THERMAL SMOOTHING OF ROUGH SURFACES IN VACUO 


Georg Wahl 
Heidelberg 

Rough metal surfaces are not stable at high temperatures: 
their structure alters over time . Lengthy baking generally re- 
sults in a smoothing of the original roughness. The fact that 
grooving at grain boundaries may occasionally result in new 
irregularities is not discussed here (we refer to the biblio- 
graphy [1] ) • The mechanism of surface alteration is presently 
understood by the concepts of W. W. Mullins [1-3] . In this paper, 
the equations for the smoothing are derived and the most important 
constants are given. 

1. The Behavior of Surfaces in Thermodynamic Equilibrium 

Before describing the time behavior of unstable surfaces, we 
shall briefly go into the thermodynamics of a surface. The most 
important thermodynamic quantity is the surface tension [4] : 


y = SF/5A ( 1 ) 

Here, 9F indicates the work that must be done to increase the area 
A by 9A at constant pressure, constant temperature and constant 
volume . 


* 

Numbers in the margin indicate pagination in the original text. 
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The surface tension of solids is between 1000 and 3000 erg cm , 
i.e., greater than that of liquids by a factor of 10. It is only 
poorly understood that this is very sensitive to impurities and 
also depends on the temperature and the orientation. Figure 1 
shows the orientation function of tungsten for the [Oil] -zone 
[5]. This is very slight, and shall therefore be disregarded 
hereafter. 


w 



[Hem bezagene Oberflachenspannung y/y m ( a 


Fig. 1. Surface tension of W for the 
[Oil] -zone [5] : y-qq “ surface tension 

in the direction 110, used as the 
abscissa (other directions 100, 111, 112 
and 332 are shown as radii) ; y - surface 
tension in any given direction. 

Key: a - adjusted surface tension y/y 


The equilibrium form of a solid at constant volume is given 
by the condition [1] : 


$ J y dA = 0 


( 2 ) 


If equation (2) does not apply, the surface is not in equilibrium: 
it changes over time. 
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2. Characteristic Roughness Functions 


There are several possibilities of mathematical description 
of roughness: 

1. The roughness function z = f(x,y) gives the position of 
the surface (in terms of the height z above the reference plane) 
at any given point of the surface (with coordinates x and y) , 
relative to a fixed reference plane. This function is very suit- 
able for a graphic description. 

2. The roughness can also be characterized by the vapor 
pressure, which depends on, the surface curvature. We have [1]: 


RT In (p/po) = y Vm K ( 3 ) 

with R as the general gas constant, T as the Kelvin temperature, 
V M as the molar volume of metal, or the volume specific to the 
quantity of matter, Pq as the vapor pressure of the regular sur- 
face and p as the vapor pressure above the curved surface. The 
factor K gives the average curvature multiplied by -2 [6] , which 
for slight inclinations can be expressed by the Laplace operator, 
applied to the roughness function f: 

* = (4) 

Thus, equation (3) is simplified: 

RT In (p/po) = — y Km Af (5) 


3. Instead of the curvature-dependent vapor pressure, we can 
use the chemical potential y corresponding to the particular vapor 
pressure: 

, (I — po + RT]n (p/p 0 ) 


( 6 ) 
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where y^ is the chemical potential of the equilibrium vapor pres- 
sure p Q above a regular surface. Between y as per equation (6) 
and Af, according to equations (3) and (4), the relation obtains: 

(7) 

/< = W - y Af ' ' 


3 . The Smoothing Equations 


Three mechanisms are responsible for alteration of surface 
structures: 1) the mechanism of evaporation, 2) volume diffusion, 

3) surface diffusion. In order to derive the smoothing equation 
for the evaporation mechanism, we proceed from equation (3) and 
calculate the curvature-dependent evaporation current density i: • 


V2nMRT 


Po oV *"■ ■*/* T 


( 8 ) 


where M is the molar mass (assuming an evaporation coefficient 
a = 1) . The evaporation current density i produces an eventual 
change in the roughness function f: 


where t is the time. 


For slight curvatures, K can be expressed by equation (4) 
and equation (8) can be developed to the second term, which gives: 


l 6f Po /. I'm v A t \ 

V K dt ~ \/2nMRT\ RT ) 


( 10 ) 


The first term on the right of (10) is equal to the evaporation 
current density over a regular surface and produces only a parallel 
displacement of the surface, but no smoothing. The smoothing 
proper is achieved only by the second term. If we ignore the 
parallel displacement, which has no significance to the smoothing. 
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we may describe the smoothing by the equation: 


3f 
8 1 


Po I'm y 
fonMRT RT 


Af 


which was first derived by W. W. Mullins [2] . 


( 11 ) 


For a one-dimensional, sinusoidal roughness, described by: 

2 = zoBin(fcz) (12) 

where z^ is the amplitude of z and k is a reciprocal length, the 
time function of the smoothness is given by: 

z = z 0 e-< , «“‘sin(fcxF (13) 


with the evaporation parameter: 


PverO — 


Po 




)/2nMRT RT 


Jfc* 


(14) 


The smoothness equations for the volume diffusion mechanism 
have also been derived by W. W. Mullins [1,3]. He proceeded on 
the assumption that the concentration of voids C T at the metal 

J-J 

surface is dependent on the curvature and given by: 

C L = C 0 -C 0 Vti~j,K (15) 


C can be either larger or smaller than the concentration of voids 

1_J 

Cq on a regular surface. This produces a diffusion stream of voids 
between the surface and the interior and a corresponding mass trans- 
port, resulting in smoothing'*'. For a one-dimensional roughness, 


'*'For a further derivation of the integrodif ferential equation 
describing the smoothing, cf. the bib, liography [1,3]. 
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given by equation (12) , the time function can be calculated by: 

z = to ein(fcz) ( 16 ) 

with the volume diffusion parameter: 

< 17 > 

in which D , is the self-diffusion coefficient for volume diffu- 
vol 

sion. 


To derive the smoothness equations by the surface diffusion 
mechanism, following Mullins [1] we start with equations (6) and 
(7) and write the surface diffusion flow i by the Nernst-Einstein 
equation : 

, _ _ w « Da 

8 RT grad * (18) 


In (18) , n g is the surface concentration of diffusing particles 
and is an empirical diffusion constant. Equation (18) can also 
be used for volume diffusion in a three-dimensional expansion [7]. 


Both empirical quantities n g and are combined by writing: 


No D t 

* 8 liy - grad fi 


(19) 


with D = (n /N n )D, . The quantity D is known as the surface 

diffusion constant; Nq is the concentration of lattice particles 

at the surface. This concentration can be estimated from V.. and 

M 

the Loschmidt number N : 

J_l 


l^L Vi 


( 20 ) 


Using the continuity equation: 
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( 21 ) 


^ll + divi * = 0 


and taking account of equations (19) and (7), we obtain the 
smoothing equation : 


at , r . u ”o”» 

0( + l * Y RT 


AAf = 0 


( 22 ) 


The symbol AA designates the double Laplace operator. 

For a one-dimensional sinusoidal roughness, the eventual 
subsidence, according to equation (22) can be described by: 


z = zo e-^o'sin (kx) 


(23) 


with as the surface diffusion parameter 


_n 1 NoD L 

P ob — RT K 


(24) 


4. Comparison of the Smoothing Mechanisms 


Figure 2 shows, using the example of smoothing of a one- 
dimensional sinusoidal roughness on a tungsten surface, the 
temperature/wavelength regions in which the individual mechanisms 
prevail. The lines bordering the regions are given by: 


and 


Pveifl/^vol — 1 


(25a) 


Pvollfio* ~~ t 


(25b) 


The constants D vq1 and D g required for the calculation (Fig. 3,4), 
as well as the vapor pressures, have been taken from [8] . Figure 2 
reveals that, for low temperatures and small wavelength of roughness, 
the surface diffusion mechanism prevails. 
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Fig. 2. Regions in which surface diffu- 
sion, volume diffusion and evaporation 
are the prevailing smoothing mechanisms 
on a tungsten surface: A-C, region of 

surface diffusion, volume diffusion and 
evaporation; a and b, boundary curves 
as per equations (25a) and (25b); X, 
wavelength of the roughness, T, Kelvin 
temperature; T g , melting point. Key: 

1 - adjusted reciprocal temperature. 



Fig. 3. Surface diffusion constants D 
for cubic body centered metals: T, 

Kelvin temperature; Tg , melting point; 
a, curve of equation (26); other values 
of Dg taken from a table in [9] and from 
[10], specifically: b for W after B. C. 

Allen (1966), c for. Mo after B. C. Allen 
(1969), d for Mo with 0.003% C after 
B. C. Allen (1966), e for W after J. P. 
Barbour (1960) , f for Mo after R. W. 
Strayer (1965), g for W after P. Bettler 
(1965). Key: 1 - adjusted reciprocal 

temperature . 



Fig. 4. Volume diffusion constants D vq ^ 

for cubic body centered metals: T and 

T s as in Fig. 3; the values for D vq1 for 

the metals indicated on the curves taken 
from [11] (a is the curve of equation 
(27). The numbers on the lines for Li, 

Mo and Nb designate the ?jy^r~ regions for 

which measurements are available for these 
lines according to [11] . For better 
clarity, these lines are not fully shown. 
Since the trend of the lines for Li and 
Na is very similar, only the Li line is 
shown and Na is indicated ' in parentheses 
at the Li line. Key: 1 - adjusted re- 

ciprocal temperature; 2 - 1-1.35. 



Fig. 5. Surface diffusion constants Dg for 
cubic face centered metals: T and T$ as in 

Fig. 3; a, curve of equation (28); the values 
for D_ taken from a table in [9] , specifically 
b for Cu after N. A. Gjostein (1961) , c for Ag 
after G. Rhead (1965), d for Cu after F. J. 
Bradshaw (1964) , e for Ni (two measurement 
series) after P. S. Maija (1967). Key: 1 - 

adjusted reciprocal temperature. 





Fig. 6. Volume diffusion constants D 

vol 

for cubic face centered metals: T and 

T 0 as in Fig. 3; the values for D , for 
S vol 

the metals indicated on the curves taken 

from [11]; a, curve of equation (29). 

Key: 1 - adjusted reciprocal temperature. 


5. Empirical Relations for Surface and Volume Diffusion 


A simple empirical rule obtains for both surface and volume 
diffusion, which can be used to determine the diffusion coeffi- 
cients of many metals. In Fig. 3 through 6, the measured values 
of the surface diffusion coefficients and the volume diffusion 
coefficients are plotted against the reciprocal temperature T, 
adjusted to the melting point T g . It can be seen that the diffu- 
sion constants may be approximately calculated by the simple 
empirical relations. For cubic body centered metals we have: 


x>. 

cm 8 b -1 


3,6 e -10 Wr 


(26) 


■Pypl 
cm* s -1 


1,0 e-VTtlT 


(27) 
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and for cubic face centered metals: 


El - 20 e -1 *- 4 Tt ^ T 

cm 2 B - " 1 


( 28 ) 
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~r =0,45e-i85’*/7 

cm* e _1 


( 29 ) 


It was N. A. Gjostein [9] who first selected such representation 
of the surface diffusion constants. 
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